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ABSTRACT 

The spectral shape of radiation due to Inverse Compton Scattering is ana- 
lyzed, in the Thomson and the Klein-Nishina regime, for electron distributions 
with exponential cut-off. We derive analytical, asymptotic expressions for the 
spectrum close to the maximum cut-off region. We consider monoenergetic, 
Planckian and Synchrotron photons as target photon fields. These approxima- 
tions provide a direct link between the distribution of parent electrons and the 
up-scattered spectrum at the cut-off region. 
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Introduction 



The interaction of relativistic electrons with low energy radiation through Inverse Comp- 
ton Scattering (ICS) provides one of the principal 7-ray production processes in astrophysics. 
In a variety of astrophysical environments, from very compact objects like pulsars and Active 
Galactic Nuclei (AGN) to extended sources like supernova remnants and clusters of galaxies, 
low energy photons are effectively boosted to high energies through this mechanism. 



The basic features of the ICS have been analyzed by I Jonesl (119681 ). iBlumenthal fc Gould 



fliaaiq ), iNarrinen & Putanenl (11993f ). iBrunettil (120001 ) and ISazonov & Sunvaevl (I2000h . The 



1970). The case of anisotropic electr o ns and /or photo ns ha s been studied by I Aharonian fc Atoyanl 



impact of the Klein-N ishina effect on the formation of the energy distribution of electrons 
was first realized by iBlumenthall (ll97lf ). Its importance in astrophysics has been dis- 
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Generally, the energy spectrum as well as the effects related to ICS are numerically 
calculated using the exact expression for the Compton cross section. On the other hand, 
compact, analytical approximations can serve as useful tools for a deeper understanding 
of the properties of Compton radiation and the implications of the complex numerical cal- 
culations. In particular, inferring the energy distribution of the parent particles from the 
observed spectrum is a much more efficient procedure when analytical approximations are 
available. For example, a power-law distribution of electrons normally results in power-law 
photon spectra. If the observed photon index is Y (in a dN 1 /dE 1 oc E~ T representation), 
then we can obtain the power-law index of the electron distribution dN e /dE e oc E~ Te from 
the relation r e = 2T — 1 for the Thomson regime and r e V — 1 if the scattering occurs in 
the Klein-Nishina regime. 

This however only applies to the energy interval far from the cut-off, the " main" part of 
the electron distribution. At the highest (and lowest) energies, there should be a break/cut- 
off in the electron distribution and in fact, the corresponding break at the radiation spectrum 
contains a lot of interesting information on the parent electrons. In particular, the peaks in 
the Spectral Energy Distribution (SED) appear at this energy range in the majority of cases, 
indicating that the source luminosity is mostly radiated at the maximum cut-off. Moreover, 
as the main, power-law part of the distribution, the shape of the cut-off carries as well 
important implication for the acceleration of the particles and in general the mechanisms 
acting in the source. Although the shape of the spectrum close to the highest energy cut-off 
is critical, this topic has not yet been adequately addressed. In this paper we examine the 
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shape of the Compton spectrum close to the maximum cut-off and we derive convenient 
analytical formulas that allow to approximate the radiated flux in this specific energy range. 

In general, the shape of the electron distribution around the cut-off can be expressed 
as an exponential, exp[— (E e /E c )^]. This allows us to describe a quite broad range of dis- 
tributions, even very sharp, abrupt, step-function like cut-offs for /5 > 1. Apart from 
the convenience of such a mathematical description, exponential cut-offs naturally arise in 
theoretical considerations. For example, in diffusive shock acceleration, power-law particle 
distributions with expone n tial cu t-off are formed when (synchrotron) energy losses are taken 
into account IWebb et al.l (119841 ) and the cut-off index is very important for investigating 



the acceleration mechanism. Analytic solutions show that in the case of Bohm diffusion 
a simple exponential cut-off exp [— E e /E c ] arises, whereas the index approaches — 2 if 
E e oc El type energy losses are taken into account, e.g. synchrotron or Thomson losses (see 
Zirakashvili fc Aharonianl 120071 ). 



In stochastic acceleration scenarios, where pile-up particle distributions are formed when 
acceleration is balanced by synchrotron type losses, th e shape of the electron cut-off is 



direc tly related to the assumed turbulent wave spectrum (I S chli ckeiserl 1 1 9 8 5l lAharonian et al. 
19861 ). e.g. (3 = 5/3 for Kolmogorov, (3 = 3/2 for Kraicman like or /3 = 2 for the hard 



sphere approximation. Of course, if more complicated energy losses dominate, like in the 
case of Klein-Nishina losses in radiation-dominated environments, more complex shapes for 
the electron distribution c ut-off may be expected in both stochastic and diffusive shock 
acceleration scenarios (e.g. IStawarz fc Petrosianll2008L IVannoni et al.ll2009l ) 



Nevertheless, it seems reasonable to consider particle distributions that exhibit expo 
nential cut -off in a general form, for investig ating and modeling the radiated spectra. In iFritz 
( 119891 ) and IZirakashvili fc Aharonianl ( 120071 ). the shape of the synchrotron spectrum close to 
the high energy cut-off has been discussed. They found that when the electron distribution 
possesses an exponential cut-off of index (3, exp [—(Ee/Ec) 13 ], then the radiated synchrotron 
spectrum exhibits a smoother cut-off, of index (3/ ((3 + 2). Apart from the practical impor- 
tance, this analytic result demonstrates that a 5-function approximation for the synchrotron 
radiation emissivity does not give the correct result. 

Here we examine the corresponding Compton spectrum, in the Thomson and Klein- 
Nishina regimes, considering different target photon fields so that both synchrotron self- 
Compton (SSC) and external Compton (EC) scattering can be addressed. We derive ana- 
lytically the asymptotic behavior of the up-scattered photon distribution close to the cut-off 
region. We consider a general electron energy distribution of the form 



dNe 
dE K 



F e (E e 
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where E e = jmc 2 is the electron energy and E c = r y c mc 2 is the cut-off energy. This presen- 
tation allows us to consider either a power-law distribution (for a = —\oc\) with exponential 
cut-off or a relativistic Maxwell-like distribution (for a = 2), that may be formed in stochas- 
tic acceleration scenarios. We consider monochromatic, Planckian and synchrotron photons 
as target photon fields. The resulting IC spectral shape is discussed for the Thomson and the 
Klein-Nishina regime and we show that is not always identical to the synchrotron spectrum, 
as is often silently assumed. Finally, we discuss basic features and physical properties of the 
radiated spectrum. 



2. Compton spectrum for monochromatic photons 



In this section we calculate the asymptotic behavior that the Compton spectrum ex- 
hibits close to the maximum cut-off, when monochromatic photons (with isotropic angular 
distribution) are up-scattered. Let as consider the general function of eq. (JTJ) , that describes 
the differential number of electrons. Electrons are considered to be isotropically and homo- 
geneously distribu ted in space. Then the spe ctrum of photons generated per unit time due 
to ICS is (see e.g. iBlumenthal fc Gouldlll~970f ) 



dN^/dE^ 



W/(£ e ,e 7 ,£ 7 )F e (£ e ) n p h{£f)dE e de 



7) 



(2) 
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V 



4e 7 E e 
m 2 c 4 



Q 



(3) 
(4) 

describes the total scattering probability, taking 
into account Klein-Nishina effects. The parameter rj in eq. (jSJ) defines the domain of the 
scattering. For rj << 1 the Thomson regime applies whereas for rj ^> 1 we are in the Klein 
Nishina regime. In the case of monochromatic photons, the number density is n p / l (e 7 ) = 
no^(e 7 — £q) and we will set no equal to 1 for this case. From now on we set mc 2 = 1 
throughout the calculations for simplicity, apart from the formulas at which the results are 
demonstrated. The case of a < will be referred to as power-law distribution, whereas 
a = 2 will be referred to as Maxwellian distribution. 



Here the function W(E e , e 7 , £7 7 ) in eq. 
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2.1. Thomson Regime 

In the limiting case oi4eoE c <C 1, when all the scatterings occur in the Thomson regime, 
the photons take a small fraction of the electron energy. Thus, from the previous relation, 
it follows that £? 7 <1/ (4eo) an d the lower limit of the integration becomes 



E e > E emin « y/E^/(4£ ) > E 1 . (6) 
Therefore, in this case rjq = E^/E e «C 1 and eq. (jl]) can be written as 

W(E e , e , E 7 ) = [2q In q + (1 - q) (1 + 2g)] , (7) 

-C / 7 max 

and in this approximation q = E^/E imax , where -E 7max = ^^E 2 . Using the electron distri- 
bution of eq. (0Q) and changing the integration variable from E e to q, the integral of eq. (j2j) 
for the Compton spectrum becomes 

f 1 2nr 2 c £ a +i 

dN,/dE,= J -^-e 2 £ 7 2 e ^ 2 g-— [2g lng + (1 - g) (1 + 2q)\ dq, (8) 

where 



We are interested in the behavior of the spectrum near the exponential cut-off, i.e. E^ » 
or £ >> 1- Then the integrant is dominated by values of q very close to unity and in 
order to perform the above integration it is convenient to change again variables to r = q~^l 2 
so that 

dN,/dE,= / — 2 £ 7 2 /(r)e-^dr, (10) 



where 

a — 3 a— 1 a — 3 a — 5 

/5 T 

The function /(r) is also dominated by values of r close to unity. Thus, we expand /(r) 
in series around r = 1. By keeping terms up to first order, the resulting spectrum reads 



dN y 



dE^ 



tK + l Ct — 1 



-e 



The above expression gives the asymptotic behavior of the Thomson spectrum with expo- 
nential accuracy at the cut-off region. When monochromatic photons are up-scattered by 
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electrons with index /3, the radiated flux exhibits a smoother cut-off, of index (3/2, i.e. 

T 



dE~ 



oc exp 



Ery[mc / 
4e n E? 



2\2\ P/2 



(13) 



As expected, the cut-off in the photon spectrum, 4e E 2 /(mc 2 ) 2 , corresponds to the maximum 
photon energy that an electron of energy E c can radiate in the Thomson regime. We also 
note that in this case, an abrupt cut-off (/3 — > oo) of the electron energy distribution would 
correspond to an abrupt cut-off of the photon spectrum. In figures [T] and [2] the analytic 
formula of eq. ( TT21) and the full, numerical spectrum are plotted, for a Maxwellian and 
a power-law distribution respectively. Asymptotics are better for a = 2. For the power- 
law distribution the numerical and analytical solution converge for very large values of the 
parameter £. In both cases, very close to the cut-off energy, the numerical spectrum is 
smoother than the approximated one. 

Finally, for a "pure" power-law distribution without exponential cut-off (a negative and 
{3 = 0) one can integrate eq. (fit)]) directly to obtain the result 

T 

27r r 2 e cA(mc 2 f szl 4(a 2 -4a + ll) 



dE^ 



2° 



£ 7 2 



la 



3) 2 (a 2 



Qa + 5) 



if a < 3, 



(14) 



that demonstrates that in the Thomson regime the radiated spectrum follows a power-law 

a-l 

of the form dN 1 / dE^ oc E y 2 . 



2.2. Klein-Nishina Regime 



When 7] = 4EoE e 3> 1, photons take almost all the energy of the electrons in one scattering. 
Then, we may define the parameter ( = AeqE^ and the low limit of the integration in eq. 
([3]) becomes E emin E 1 . Let us change variables to x = E^jE e . As ( ^> 1, we keep only 
the leading terms 



dNry 

dE^ 



8nr 2 c 



i-i/e 



x- 



2(1 



x) 



E e [E 1 jx)dx. 



(15) 



Here we have approximated the upper limit of the integration by 1 — 1/C because the integral 
in eq. f JT5|) diverges logarithmically at x — 1. For an electron distribution with exponential 
cut-off the main contribution to the integral comes from regions of x close to unity, so that 



1 < 



2(1-3=) 



and we can neglect the first term from the expression in the brackets, 



JeZ 



Anr 2 c 



x 2 dx 



x 
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X 



(16) 
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Now the integral of eq. ffl6|) can be calculated for energies close to the cut-off region, 
E y ^> E c , leading to a spectrum of the form 

KN 



tZJVL 



£()Ey 



In 



mc 



.2\2 



In /? - P ln( 



-7 



F e (Ey), for £ 7 >£ c , 



(17) 

where 7 = 0.5772 is the Euler' s constant. In the Klein-Nishina regime where electrons loose 
almost all their energy in each scattering, the Compton Spectrum practically reflects the 
behavior of the electron distribution. Thus, in this case, the exponential cut-off maintains 
the index /3 and is steeper than in the Thomson case 

KN 



dE~, 



oc exp 



E, 



P 



'11 



As before, the photon cut-off energy corresponds to the maximum photon energy that elec- 
trons of energy E c radiate in the Klein-Nishina regime. Moreover, an abrupt electron dis- 
tribution cut-off would result in an abrupt photon spectrum cut-off. The asymptotics of eq. 
( 117]) are plotted in figures [3] and IU as can be seen, they provide a good approximation just 
after the peak of the SED. 

The fact that the exponential index becomes 0/2 and /3 in the Thomson and Klein- 
Nishina regime, respectively, indicates that using a 5-function approximation for the cross 
section provides a correct result for the calculated spectrum in these two regimes (for an 
exten ded discussion on the applicability of <5-function approximation see ICoppi fc Blandford 
19901 ). Obviously this is not true for values of r\ c close to unity (see fig. EJ), where r\ c = 4:EqE c 
refers to the electron cut-off energy. As we are interested in the highest energy part of the 
spectrum, the approximation in the Klein-Nishina regime is satisfactory even for values of 
T) c that do not significantly exceed unity. This happens because since r) c > 1 for the electron 
cut-off energy it holds for all the energies E e > E c that actually form the shape of the 
exponential cut-off. On the contrary, for the Thomson regime one needs all the radiated 
photons above the cut-off to be emitted at this regime, which indicates rather small values 
of 77 for the approximation to be good, especially for low (3 factors. 

A similar calculation can be performed for a "pure" power-law electron distribution, 
with (3 = 0. If we rewrite the second term of eq. ffl5|) in the form 



/ x- a F e (E,/x) = F e (Ej) / + F e (Ej) 

Jo 1 - x J 1 - x J 



X 



'0 x — x JO x ~ x ^0 

we can perform the above integration resulting in an emitted spectrum of the form 

KN 

■Kr 2 cA(mc ) 



dE^ 
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(19) 



(20) 
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where ty(a) is the digamma function defined as the logarithmic derivative of the T func- 
tion, \l/(a) = r'(a)/r(a). This formula shows that the emitted compton spectrum is much 
steeper in this than in the Thomson regime due to the suppression of the cross-section. The 
functional dependance of eq. (j2D]l 



dAL 



KN 



oc E 



-a-l 



(21) 



has bee n obtained in lBlumenthal fc Gouldl (jl970l ). see e.g. their eq. (2.87), and lAharonian fc Atoyan 
(1l981o ) . see their eq. (32). These formulas differ in the term related to the power- law index 
a, due to the different approach used for the calculation of the asymptotics. This difference 
is negligible. 



3. Compton spectrum for a broad photon distribution 

Once we have calculated the radiated spectrum for monochromatic photons, we can 
examine the behavior of the Compton spectrum for various photon fields. The case of 
monochromatic photons is important for understanding the scattering mechanism and a 
necessary step for further calculations. However, in nature the photon fields are usually 
broader than the monochromatic one, except if we deal with emission lines. Here we will 
consider a Planckian photon distribution which is often the case in external Compton sce- 
narios and we will examine as well the case of synchrotron photons from the same parent 
electron distribution, that are used as the target photon field in synchrotron self-Compton 
models. 



3.1. Planckian photon field 

Let as assume a Planckian distribution for the photon field so that the differential 
number density is given by 

1 ^ 

For the Thomson regime, where now we demand 4kTE e <C 1, we can use the Wien limit 
(e 7 ^> kT) at which 



1 



n*(ey) = -^T-A^' kT - (23) 



7T 



l h 3 c 3 



This is acceptable as the asymptotic behavior of the Compton spectrum at high energies is 
mostly defined by the soft photons with energy around and greater than kT. 
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In that case the integration of the spectrum eq. (1121) over the photon energies can be 
evaluated by the saddle point method of integration. After replacing e with e 7 in eq. (fl2]l . 
the integral is written as 

T 



dE~, 



BB 



dE~, 



poo 

/ g(E 7 ,e 7 )e s ^ 
Jo 
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where 



and 



g{E. 



7' e i) 



llXr 2 A 



E 7 2 



3-a 

e 7 



'7T 



] h 3 c 2 



S(E, 



7' e U 



E~ 

a^eI 



v 2 e 

0/2 

_ _L 
kT' 



(24) 



(25) 



(26) 



Here BB stands for black body. Let as define the parameter £i that is related to the cut-off 
of the up-scattered photon spectrum for the case of monochromatic photons, if we replace 



e 7 with kT 



6 = 



(27) 



4kTE 2 ' 

This parameter simply describes the outgoing photon energy normalized to the maximum 
energy which E c electrons radiate when they up-scatter photons of energy kT. One can 
use the saddle point method for the above integration because the integral at large ener- 
gies, £i 3> 1, is determined by the soft photon energy interval around the energy x which 



maximizes/minimizes the function S(E, 



7' e ~fJ 



The saddle point xq is at 



x 



2 



2 A 

(3 



P+2 



kT. 



(28) 



Then, the Thomson spectrum for a Planckian photon distribution is calculated to be dN 1 / dE^ 
g(E y , e*) exp[— S(E^, x )} \j2n / — S"(E 7 , x ) where S" the second derivative of S at the sad- 
dle point xq. After rearranging the terms, we retrieve the following expression 

T 



dE^ 



Anr 2 e A(kT) 



2 mc 



2\a+l 



BB 



2 a 7l 2 H 3 C 2 



71 



(3 + 2 



X r. 



7 e: 



-£+2(2|l )^+2 



for 



4kTE 2 

7 ^ I T\1 ' 



(29) 



Therefore, when the target photon field is a black-body, the shape of the cut-off is 
affected by the soft photon distribution and the exponential cut-off in the Thomson spectrum 
is smoother in comparison to the monochromatic photons case. The index now becomes 
0/(/3 + 2) as 



diV 7 
dE^ 



cx exp 



(3 + 2 



BB 



2 E^mc 2 ) 2 
/3 AkTE 2 



d + 2 



(30) 
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This exponential cut-off is always smooth (less than unity) and it becomes unity in the 
case of an abrupt electron distribution cut-off, as lim^^oo (3/ ((3 + 2) = 1. Interestingly, the 
Thomson spectrum for Planckian photons at high energies exhibits the same exponential cut- 
off shape as the synchrotron spectrum. This is the only case where the two components of the 
spectrum show the same behavior for arbitrary index (3. For the Maxwellian and power-law 
type distributions of electrons, Eq. [30] is presented in fig. [6] and fig. [7] respectively. 

In the Klein-Nishina regime the shape of the exponential cut-off does not depend on 
the up-scattered photon distribution but, as mentioned above, preserves the electron index 
(3. Integration of eq. f|T7|) over photon energies (after replacing e with e 7 ) requires the 
calculation of the following integrals 

r°° i (kT) 2 

/ — n(e 7 ) In (4e 7 £ 7 )de 7 = (In (AkTE^ - 0.1472) . (32) 

Jo £7 bfic 

Then, the asymptotic behavior of the up-scattered spectrum close to the cut-off follows the 
formula 



dN 1 



dE^ 



KN 



7T 2 r 2 e A(KT) 2 (mc 2 ) 2 F e (£ 7 ) 



W>c 2 £ 7 

BB ' 



ARTE E 
ln^-^-ln/3-/31n^-0.724 

{mc z ) z E c 



for £ 7 > E c 
(33) 



Thus, in the Klein-Nishina regime the cut-off is always much sharper than in the Thom- 
son regime. The spectra for Maxwellian and power-law distributions of electrons are shown 
in fig. [H]and fig. M respectively. For values of the index (3 = 1,2,3, the corresponding 
shape in Thomson regime becomes 1/3, 1/2 and 3/5 respectively, always less than unity. For 
an abrupt cut-off {(3 — > 00), the Klein-Nishina spectrum appears sharp as well, while the 
Thomson spectrum exhibits a simple exponential cut-off. 



3.2. Synchrotron photon field 

In synchrotron self-Compton models, the electrons up-scatter the photon which they 
produced via synchrotron radiation. In contrast to external Compton models, we do not have 
an analytic expression for the target photon density. However, we can use an approximation 
for the synchrotron spectrum at energies around the synchrotron cut-off, i.e. for e 7 > bE 2 , 
given that the main contribution to the scattering process at high energies comes from this 
energy range. Here b = 3qBh/4irmc(mc 2 ) 2 . 
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3.3. Synchrotron spectrum 



In the case of chaotic magnetic fields, the synchrotron emissivity of an electron with 
energy E e is described by the equation 



dN-y _ V3q 3 B g /ey 
de^ mc 2 he^ 



where 



6 7 2 



3qBh E 2 e 
Anmc (mc 2 ) 2 



(34) 



(35) 



is the "critical" energy for synchrotron em ission. The synchrotron power e 7 iV 7 /<ie 7 peaks 
at 0.29e s (e.g. iRybicki &: Lightmannl Il979l ) The function G(y) can be well approximated, 
with an accuracy bett er than 0.2% over the entire range of variable y, by the formula 
( lAharonian et al.l 120101 ) 



G(y) = G{y)e~y 



I.8O82/ 1 / 3 1 + 2.21i/ 2 / 3 + 0.347?/ 4 / 3 
y/1 + 3Ay 2 / 3 1 + l-353y 2 / 3 + 0.217t/ 4 / 3 



(36) 



For large y 1, the function G(y) is approximately G(y) ~ a/vt/2. Let us define the 
parameter £2 that defines the synchrotron photon energy normalized to the energy e s (eq. 
351) for the electrons with energy E c 



62 = e 7 /e 



bE 2 



(37) 



In order to calculate the emitted synchrotron spectrum, we need to integrate eq. ( 134)) over 
the electron distribution given by eq. ([1]). By changing variables from E e to y = e 7 /e s , we 
find 



SYN 



^F e (E e )dE e = ^ q l B u A C-l)^ 
o cte 7 2mc I he 1 b 



y-^- 1 G(y)e- y -^ )m dy. (38) 



For £2 > 1 the integration over y reveals a saddle point at 



/2f 2 



2/0 = 2 V /3 





/3+2 



(39) 



so that finally the emitted synchrotron spectrum can be expressed as 

V3q 3 B 

SYN 



e^dNj 



de~ 



(40) 
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This formula indicates the shape of the exponential cut-off for synchrotron radiation 



e-ydN 7 



de~ 



oc exp 



SYN 



+ 2 (2e,\^ +2) 



(41) 



We note that the 0/( + 2 ) index for the synchrotro n expo nential cut-off has already been 
found in iFritzl (119891 ) and IZirakashvili fc Aharonianl (120071 ) . With the above calculations 
we can moreover estimate critical aspects of the emitted spectrum, and in particular the 
cut-off energy of the emitted synchrotron spectrum (see section 4). The above equations 
correspond to optically thin synchrotron sources when the synchrotron-self absorption can 
be ignored. This could be the case of even very compact and highly magnetized sources, 
as long as the synchrotron cut-off appears at optical and higher frequencies. One should 
also ment ion that the synchrotron spectrum is sensitive to inhomogeneities of the magnetic 



field (e.g. iKatz-Stone fc Rudnicklll994l . lEilek et al.lll997l ). Obviously, the fluctuations of the 



magnetic field should have an impact on the synchrotron spectrum, namely they will mak e 
the cut-off smoother and shifted towards higher frequencies (see e.g. lEilek &: Arendtlll996l ). 
In this regard, the ICS 7-ray spectrum is free of uncertainties related to the magnetic field 
distribution, except for realization of the SSC scenario in the Thomson limit. 

Last, we should mention that using a 5- function approximation for the synchrotron 
emissivity would result in an exponential cut-off of index 0/2, whereas the correct value is 
0/(2 + 0). The synchrotron asymptotics are plotted in figures [TU1 and [TTj for Maxwellian 
and power-law electron distributions, respectively. For comparison, the Thomson spectrum 
for Planckian photons is shown as well. 



3.4. SSC spectrum 

Now we can integrate the Compton spectrum for monochromatic photons over the 
photon energies for the synchrotron distribution. The differential photon number density 
(for a spherical source) is 

, , 3 R dN^ , , . 

n '^> = «?7 &7< (42) 

where R is the source size. Let us first perform the calculations for the Thomson regime 
described by eq. ( fT2|) with eq — )■ e 7 . In that case the exponential factor in the integrant is 
exp [— S], where S = £ + (2 — 0)y o /0. The function S has extremum at the saddle point 

0(2E 1 \^ 

z o = 77 Tin e o> ( 43 ) 



2 \0E O 
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where Eq = AbE^ is the maximum photon energy that results from electrons with energy 
E c when they up-scatter synchrotron photons with energy e . The analogy with the case of 
Planckian soft photons is direct. Replacing the synchrotron characteristic energy bE 2 c with 
kT results in the same saddle point. 



Here the second derivative of the exponential argument at the saddle point has a simple 



form 



S"(z ) 



(3E 1 
2E,eV 



so that the integration over the synchrotron number density gives 

T 



dE^ 



ssc 



AncR 2 



27r' z r z cA z (mc z ) a+1 V3q 3 B ( E^ 

T 



2V/3(/? + 2)/> 



G(y ) _ « 3 _i±± y ~ 



(44) 



(45) 



where y is calculated at the saddle point 



~ (3 
Vo= 2 



2E 7 
Wo 



(3+4 



Eq 



AbEt 



[mc 



2\2 



(46) 



In the case of SSC radiation, the electron distribution up-scatters the synchrotron pho- 
ton distribution with an exponential cut-off exp[— (e 7 /e c )?+ 5 ]. The corresponding Compton 
flux at high energies exhibits a cut-off index (3/(j3 + 4) which is smoother than the seed 
synchrotron distribution, 



dE^ 



oc exp 



/3 + 4 (2E, 



ssc 



\(3E C 



p 

rl + 4 



(47) 



For (3 = 1, 2, 3, the corresponding values for the Thomson exponential index is 1/5, 1/3, 3/7, 
significantly less than the electron distribution index and different than the synchrotron case. 
If /3 — > oo, then the SSC spectrum shows a simple cut-off, like in the case of up-scattering 
Planckian photons. The asympotic formula of eq. (l4"5j) is plotted in figures [12] and [13] for 
Maxwellian and power-law electrons. 

In the Klein-Nishina regime, the integration over the synchrotron photon density can 
not be performed analytically for arbitrary values of the indexes a and (3. In this regime 
however, the soft photon field does not play an important role in the shape of the up-scattered 
spectrum close to the maximum cut-off. Thus, eq. (117)) for monochromatic photons offers 
a rather good description of the asymptotic behavior of the SSC spectrum at Klein-Nishina 
regime. See e.g. the figures [TH and [15] for Maxwellian and power-law electrons respectively. 
The analytic formula describing the asymptotes is eq. ffTTl) normalized to the numerical 
solution, where the soft photon energy e 7 has been replaced by e s = bE 2 c . 
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4. Discussion 

The main results of this paper are summarized in table [TJ In the Klein-Nishina regime 
the up-scattered Compton spectrum exhibits the same exponential cut-off index (3c, as the 
electron distribution index (3, and does not depend strongly on the target photon field. This 
implies that from the 7-ray spectrum we practically "observe" the electron cut-off shape. In 
particular, an abrupt electron distribution cut-off ((3 — > 00) would result in an abrupt cut-off 
for the photon spectrum (J3 — > 00). The case of up-scattering monochromatic photons is 
shown in figures |3] and |U The case of Planckian photons is plotted in figures |8] and [9j whereas 
figures E] and [13 correspond to the SSC spectrum. In all of these cases, the asymptotics are 
rather good just after the peak of the SED. 

On the contrary, in the Thomson regime the up-scattered photon exponential shape 
is always smoother than the electron distribution cut-off shape. For monochromatic target 
photons this is flc — (3/2, as shown in figures [T] and |2J In this case, the shape of the cut-off in 
both Thomson and Klein-Nishina regime (eqs. [T51 and ITS]) shows that using a 5-function for 
the Compton emissivity provides the correct result. Here the asymptotic analytic expression 
of eq. ( !T2|) approaches better the numerical solution for a Maxwellian electron distribution. 

Interestingly, for Planckian photons we find a different relation between (3c and (3, 
which is the same as for synchrotron radiation, (3c = (3/ {(3 + 2). As can be seen from figures 
|6] and the approximation is very good. Finally, for synchrotron photons it holds that 
Pc — + 4). As in the case of monochromatic photons, the SSC asymptotics in the 
Thomson regime are better for Maxwellian electrons. 

In general, although in the Klein-Nishina regime the Compton spectrum preserves the 
electron distribution index (3c = (3, in the Thomson regime the up-scatter photon cut-off 
index is always smaller than the electron distribution cut-off, (3c < (3. The only exception is 
for monochromatic photons up-scattered by electrons with (3 — > 00. In this case the Compton 
spectrum should exhibit as well an abrupt cut-off. For Planckian and synchrotron photons 
(3 — > 00 for electrons means a simple exponential cut-off for the Compton SED, (3c = 1. 

The Thomson spectrum for Planckian photons exhibits a cut-off index (3s = (3/ ((3 + 2), 
same as the synchrotron spectrum. Thus, if the SED of the observed object is consid- 
ered to consists of these two components, synchrotron radiation for the low energies and 
up-scattering of Planckian photons in the Thomson regime for high energies, then the expo- 
nential part of the two "bumps" is very similar (see fig. [TUl and ITTj) . If, on the other hand, 
the up-scattering occurs in Klein-Nishina regime in the energy band of the cut-off, then this 
is not true. While (3$ = (3/((3 + 2) for the synchrotron component cut-off, (3c = (3 for the 
high energy component. Even if we consider an abrupt cut-off for the electron distribution, 
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the two bumps would be different ((3s = 1 and /3c — > oo, respectively). For an SSC model, 
the two components do not show the same exponential cut-off shape neither in the Thomson 
(where we get (3s = (3/((3 + 2) and (3c = (3/((3 + 4) for low and high energies respectively), 
nor in the Klein-Nishina regime ((3s — /?/(/? + 2) and (3c = (3). Only if the electron distri- 
bution has an abrupt cut-off, then in the Thomson regime we can get (3c = (3s = 1, while in 
the Klein-Nishina regime (3c — > oo. Thus, the two components of the SED do not show in 
general the same shape at the cut-off region. 

Apart from the shape of the up-scattered, photon spectrum at the cut-off region, another 
interesting point concerns the cut-off energy itself, E ltCUt . We have shown that in the Klein- 
Nishina regime this remains always at E yiCUt = E c , same as the electron distribution cut-off 
energy, independently of the soft photon field. On the contrary, at the Thomson regime, the 
cut-off photon energy depends on the target photon field. For the monochromatic target 
photons the value of E... is obvious; it is simply equal to the cutoff energy of electrons. 
For monochromatic photons, the resulting value of E lfiUt is rather obvious, simply because 
there is a maximum up-scattered photon energy for the fixed electron energy. Consequently, 
the cut-off energy of the IC spectrum must be equal to the maximum up-scattered photon 
energy by electrons of energy E c , as it follows from (eq. US]) . 

E ^ cut ~ J^f (48) 

On the other hand, in the case of a broad distribution of target photons, the relation be- 
tween E c and E 1;CUt depends on the seed photon spectrum and the index (3 of the electron 
distribution. For the Planckian photon distribution, we find (see eq. 130]) 

^ cut (mc 2 y\(3 + 2j 2" { } 

In analogy to the monoenergetic photons case, the maximum energy at which electrons of 
energy E c can radiate, when they up-scatter photons of energy kT, is 4E 2 KT/(mc 2 ) 2 . In 
respect to this, the Thomson spectrum cut-off energy is smaller by a factor of ((3/2)[2/((3 + 
2)](/ 3 + 2 )/ 2 . This factor is not negligible especially for small (3, e.g. it takes values of ~ 0.15, 
0.25, 0.33 for (3 = 1, 2, 3 respectively (almost one order of magnitude for a simple exponential 
cut-off). As expected, it tends to unity for (3 — > 00. It does not depend however on index a 
of the electron distribution. 

For synchrotron radiation, from eq. fT4T]) . we find the cut-off energy 

eS^(^Vf, (50) 
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which reveals exactly the same factor as in the case of Thomson spectrum for Planckian 
photons, but now in respect to the characteristic energy bE 2 c . 



Finally, for the case of SSC, eq. ( 117]) gives 

4bEt ( 2 \ t T B 4E 2 ^ SYN 



rpSSC 



(mc 



2\2 



(3 + 4 



(mc 



2\2 



(3+4: 



d+4 

5 



(3 + 2 



(51) 



If we compare the cut-off energy with e.g. 4bE 2 ut e^J/(mc 2 ) 2 then the factor related to 
the index (3 takes values of ~ 0.035,0.15,0.25 for (3 = 1,2,3, slightly less than in the case 
of a Planckian target photon field. This is due to the fact that the cut-off shape is now 
smoother. These analytic results are useful when modeling the observed spectrum, as one 
may infer the electron distribution cut-off energy from the photon spectrum cut-off energy 
with only the uncertainty introduced by a (possible) Doppler boosting. The position and the 
amplitude of the synchrotron and IC peaks in the SED contain very important information 
about physical parameters of non thermal sources, like the strength of the average magnetic 
field and the energy density of relativistic electrons. The shape of the SED, especially in the 
region of the cutoffs of the synchrotron and IC components of radiation, provide additional, 
more detailed information about the distributions of electrons and magnetic fields. For 
example, the spectral cutoff in the IC component formed in the the Klein-Nishina regime 
provides direct, model-independent information about the energy spectrum of highest energy 
electrons. This is a critical issue for understanding of particle acceleration mechanisms. 
Furthermore, combined with the shape of the synchrotron cut-off, it can allow us to extract 
information about the distribution of the magnetic field. This can be demonstrated by the 
following simple example. Let as assume that we have observed a smooth synchrotron cut-off 
which can be interpreted as the result of an electron distribution with an exponential index 
e.g. 8 p» 1. This hypothesis can be checked by the shape of the cut-off of the IC component. 
If the latter is formed in the Klein-Nishina regime, and exhibits a sharp cut-off behavior 
indicating to the electron distribution with 8 > 1, then one should attribute the smoothness 
of the synchrotron cut-off to magnetic field inhomogeneities rather than to the actual shape 
of the electron distribution. 



Table 1: exponential cut-off index for Compton spectrum 
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5. Summary 

In this paper we have examined the asymptotic behavior of the Compton spectrum 
close to the maximum cut-off. We assumed that the electron distribution follows the general 
formula Ef exp[— (E e /E c )^} so that our analysis may account for a relativistic Maxwellian- 
type distribution, as well as for a power-law distribution with exponential cut-off. The 
exponential cut-off of the electron energy spectrum results in an exponential cut-off in the 
Compton spectrum, of the form exp[— (£' 7 /£' 7)Cttt )^ c ], with j3 c and E 1>cut the corresponding 
cut-off index and energy respectively. We show that in the Klein-Nishina regime, the cut-off 
index remains unchanged, f3 c — p. The shape of the up-scattered spectrum close to the 
maximum cut-off basically "reflects" the electron distribution and does not depend strongly 
on the target photon field. The cut-off energy also correspond to the electron distribution 
cut-off energy, E 1>cut = E c . 

In the Thomson regime, the resulting spectrum close to the cut-off is very different. 
First of all it strongly depends on the up-scattered photon field. Monoenergetic photons 
lead to a cut-off of index f3 c = (3/2, whereas Planckian photons result in f3 c = /3/(/3 + 2). 
When the up-scattered photon field is the synchrotron photon field, as in SSC models, then 
the cut-off appears extremely smooth, with an index f3 c — f3/(/3 + 4). In contrast to the 
Klein-Nishina regime, the Thomson spectrum cut-off energy E ljCUt depends not only on the 
electron distribution cut-off energy, but also on the target photon field and as well on the 
index j3. 

The obtained analytic expressions are useful for deriving the electron spectral shape at 
the cut-off region directly from the observed high energy flux. These two parameters may 
give important insight into the acceleration and radiation mechanisms acting in the source. 
Furthermore, one may use the higher energy part of the observed Compton spectrum as a 
"diagnostic tool" to distinguish between EC and SSC models, as different photon fields lead 
to different cut-off shapes. 

Acknowledgement: We would like to thank F. M. Rieger for helpful discussions and 
careful reading of the manuscript. 
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Thomson regime 
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Log (Ey4£p h E 2 cut ) 



Fig. 1. — Asymptotics in the Thomson regime for monochromatic photons. The r\ c = 
4:EqE c parameter that defines the domain of scattering for the cut-off electrons is r] c = 
0.0004. E c = 10 2 and 6q = 10 -6 both in mc 2 units. For this figure relativistic Maxwell-like 
electron distributions are used with different shapes of exponential cut-off, ft parameter. 
The exponential cut-off of the up-scattered photon spectrum posses an index of (3/2. As 
in all figures, black solid lines show the numerical spectrum, whereas colored, dashed lines 
correspond to the analytic approximation. 
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Fig. 2. — Same as figure Q] but for power-law electrons. The asymptotics approach the 
numerical solution only for E 1 ^> 4eoE 2 ut . For smaller values of the numerical spectrum 
is smoother than the approximated one. 
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Fig. 3. — Asymptotics in the Klein-Nishina regime for monochromatic photons and Maxwell- 
like electrons. The rj c parameter that defines the domain of the scattering is r] c = 400. 
E c = 10 6 and Eq = 10~ 4 both in mc 2 units. The photon spectrum exhibits a cut-off index /3, 
same as the electron distribution. The analytical approximation is in very good agreement 
with the numerical spectrum just after the peak of the SED. 




Fig. 4. — Same as figure [3] but for a power-law electron distribution. 
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Fig. 5. — Asymptotics for a Maxwell-like distribution with (3 = 2 for different values of 
the parameter 1] C that defines the domain of the scattering for electrons with energy E c . 
In the Klein-Nishina regime the approximation is good even for not very large values of 
r] c . In the Thomson regime the approximation becomes acceptable for rj ~ 0.004 <C 1. 
In the intermediate domain the exact numerical spectrum decays more sharply than the 
approximated one, indicating that while 4e _E c < 1, part of the electrons with E e > E c 
already up-scatter the soft photon in the Klein-Nishina regime. 
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Fig. 6. — Asymptotic behavior of the Thomson spectrum at the cut-off region, for Maxwell- 
like electrons up-scattering Planckian photons. Parameters used are E c = 10 2 , KT = 10~ 6 
so that r] c = AkTE c = 0.0004 1. The exponential cut-off possesses an index /3/(/3 + 2). 
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Fig. 8. — Compton spectrum in the Klein-Nishina regime for Maxwellian electron and 
Planckian photons. E c = 10 6 , Kk = 10~ 4 and rj c = AkTE c = 400 ^> 1. The spectrum 
shows a cut-off of index /3, same as the electron energy distribution. 
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Fig. 10. — Synchrotron spectrum of relativistic Maxwell-like electrons. With black crosses 
the Thomson spectrum is shown for Planckian photons and — 2, scaled to E^/5kTE^ ut 
and normalized to the synchrotron flux. At large energies, the Thomson and synchrotron 
spectra show similar shape of the cut-off and the corresponding index is 0/(/3 + 2). The 
cut-off energy for the electrons is E c = 10 2 and the magnetic field B = 1G. 
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Fig. 12. — SSC radiation at the cut-off region in Thomson regime. A Maxwellian electron 
distribution has been used, with parameters E c = 10 2 and B = 1G for the magnetic field. 
The exponential cut-off that arises is very smooth, of index /3/(/3 + 4) 
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Fig. 13. — Same as figfT^but power-law electrons. The asymptotics approach the numerical 



solution only for E~ ^> 



Very close to the photon cut-off energy, the numerical 



spectrum is smoother than the approximated one. 
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Fig. 14. — SSC radiation at the cut-off region in the Klein-Nishina regime for Maxwellian 
electrons. The parameters used are E c = 10 5 and B = 1G, so that the parameter r] c w AbE* 
that defines the domain of the scattering is ?7 C ~ 90 ^> 1. The analytical formula plotted here 
is eq. ( 112]) for monoenergetic photons with e 7 = bE 2 cuV scaled to the numerical spectrum. 
The cut-off shape has the same index /3, as the electron distribution. 
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Same as figJH] but power-law electrons. 
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